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Abstract 

Two new observations have strengthened the case for the supersymmetric nature of the Cold 
Dark Matter component in our universe: First, it was shown that new data on the nuclear abun- 
dance, B/C - and ^"Be/^Be ratios constrain the diffusion parameters in Galactic Models so strongly, 
that they lead to a clear deficiency in the production of diffuse hard gamma rays, antiprotons and 
hard positrons, if no anomalous sources or anomalous energy dependence of the diffusion coeffi- 
cients are postulated. Second, from the precise relic density measurement by WMAP the WIMP 
annihilation cross section can be determined in a model independent way. If the WIMPS are 
postulated to be the neutralinos of Supersymmetry, then only a limited region of parameter space 
matches this annihilation cross section. It is shown that the resulting positrons, antiprotons and 
gamma rays from the neutralino annihilation (mainly into bb quark pairs) provide the correct shape 
and order of magnitude for the missing fluxes in the Galactic Models. Here the shape is the most 
discriminating feature, since for relatively heavy neutralinos the positron and gamma spectra from 
neutralino annihilation are significantly harder than the spectra produced by secondary interac- 
tions. For the antiprotons the shapes of background and signal are similar, but they still provide a 
strong constraint by the requirement of not overproducing antiprotons. The fitted normalization 
factors for antiprotons, positrons and gammas from neutralino annihilation come all out to be sim- 
ilar and below 10 for a standard NRW halo profile. It should be pointed out, that it is absolutely 
non-trivial to solve all three deficiencies simultaneously with a common halo profile, because the 
mean free paths of photons, antiprotons and positrons are quite different, so the observed fluxes 
come from different regions of the galaxy. The probability of a global fit to the galactic spectra of 
diffuse gamma rays, positrons and antiprotons improves from about 10~® to 0.5, if Dark Matter, 
as predicted by Supersymmetry, is taken into account. This corresponds to about 6a evidence for 
the supersymmetric nature of Dark Matter in case of Gaussian errors. 



1 Introduction 

Cold Dark Matter (CDM) makes up 23% of the energy of the universe, as deduced from the tem- 
perature anisotropies in the Cosmic Microwave Background (CMB) in combination with data on the 
Hubble expansion and the density fluctuations in the universe The nature of the CDM is 

unknown, but one of the most popular explanation for it is the neutralino, a stable neutral particle 
predicted by Supersymmetry |S1E1- The neutralinos are spin 1/2 Majorana particles, which can an- 
nihilate into pairs of Standard Model (SM) particles. The stable decay and fragmentation products 
are neutrinos, photons, protons, antiprotons, electrons and positrons. From these, the protons and 
electrons are drown in the many matter particles in the universe, but the others may be detectable 
above the background from nuclear interactions, especially because of the much harder spectra ex- 
pected from neutralino annihilation. The background is strongly constrained by the recent, precise 
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measurements of the fluxes of all nuclei in the galaxy, especially the ratios of secondary /primary nu- 
clei and radioactive/non-radioactive isotopes. The diffusion and reacceleration parameters deduced 
from these data describe the fluxes of all matter particles in the galaxy, but they predict too few 
antiprotons, hard gammas, and positrons [71|H]. Several ad hoc proposals have been made to solve 
these deficiencies, like unphysical breaks in the diffusion coefficient or postulating "unprocessed" fresh 
components in a "Local Bubble", which is different from the rest of the galaxy [^j, or in order to ob- 
tain more high energy gammas one has either to postulate a local hard nuclear component generating 
hard tto's or a local hard electron spectrum to create hard gammas by inverse compton scattering (Hj. 

However, none of these proposals solves these deficiencies simultaneously. In this paper we consider 
the annihilation of Dark Matter particles as a source for the missing gamma rays, antiprotons and 
positrons. A statistical analysis of all three species simultaneously provides a rather stringent test, 
since the mean free path is quite different for gamma rays, antiprotons and positrons. If one can 
solve the deficiencies of all three simultaneously for a single halo density profile and a single set of 
SUSY parameters, then this can be considered as a 'smoking gun' signature for the supersymmetric 
nature of the Cold Dark Matter, at least if the statistical significance is sufficient. Fortunately, the 
recent improvement in the cosmic parameters severely constrains the background calculation, while 
the precise WMAP data on the relic density implies a precise and model-independent measurement 
of the annihilation cross section: a larger annihilation cross section would imply a lower relic density 
and vice versa. Such an annihilation cross section can only be obtained for rather restricted regions 
of the SUSY parameter space. It will be shown that these SUSY parameters combined with the best 
known Galactic Model improve the probability of the fit to the measured galactic fiuxes of diffuse 
gamma rays, antiprotons and positrons by 8 orders of magnitude, which corresponds to more than 
6a evidence for the supersymmetric nature of Dark Matter in case of Gaussian errors. 

Indirect detection of Dark Matter has been discussed much before jJOI-IEB ■ results differ from 
these previous results by performing a statistical analysis to positrons, antiprotons and gamma rays 
simultaneously and taking into account the best known propagation models. Furthermore, previous 
results never considered the shape of the gamma spectrum, which turns out to be a discriminating 
feature of the present analysis, since for the relatively heavy neutralinos the diffuse gamma spectrum 
from ttq decays is considerably harder than the spectrum from secondary interactions. Also most of 
the previous analysis were done before the WMAP data became available and before the new satellite 
data on nuclear fluxes constrained the galactic background processes. 

We do not consider neutrino detection, since for neutrinos no precise data are available. We 
restrict us to neutralino annihilation in the Minimal Supersymetric Model (MSSM) with gravity 
inspired supersymmetry breaking and radiative electroweak symmetry breaking, since this simplest 
model already gives a good description of the considered fluxes. Actually all we need from the model is 
stable neutralinos in the mass range of a few 100 GeV, which annihilate predominantly into hadronic 
final states. The cross sections turn out to be correct for values of the ratio of neutral Higgs vacuum 
expectation values around 50, as expected in SO(IO) type of models|.52_. In this case the neutralino 
annihilation into b-quarks pairs via pseudoscalar Higgs exchange dominates and the neutralino is 
predominantly photino-like, or more popular. Dark Matter is the supersymmetric partner of the 
Cosmic Microwave Background (CMB). 

In the following we first describe the model independent determination of the annihilation cross 
section from the relic density. In the next section the predictions from Supersymmetry concerning 
neutralino annihilation are discussed, while in the following section the deficiencies in the gamma rays, 
antiprotons and positrons, as predicted by the Galprop model [HSIISI] are discussed. In the last section 
the global fits are discussed. They are performed within the frame work of the DarkSusy |55j program, 
after modifying it to obtain the background and propagation model from Galprop in DarkSusy in 
order to have consistent propagation and diffusion for antiprotons and positrons from the two sources: 
nuclear interactions and neutralino annihilation. For the SUSY particle spectrum we use Suspect 
for the Higgs masses Feynhiggsfast |57| and as a cross check for the cross sections and relic density we 
use Micromegas |58) . Finally, some expectations for direct detection and indirect searches for Dark 
Matter from solar neutrinos will be given. 
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Figure 1: The running of the squark - and slepton masses starting at mo, gaugino masses starting at mi/2 
and Higgs mass parameters starting at ^ rn^ + for mo — mi/2 = 500 GeV and ta.n/3=51 (left) and for 
mo = mi/2 = 1000 GeV and tan /3=5,? (right), which are the parameters of interest for the present analysis. 



2 Annihilation Cross section Constraints from WMAP 

In the early universe all particles were produced abundantly and were in thermal equilibrium through 
annihilation and production processes. The time evolution of the number density of the particles is 
given by the Boltzmann equation, which can be written for neutralinos as: 

^ + SHn^ = _<^yy^r,l-nf), (1) 

where H is the Hubble expansion rate, is the actual number density, n'^ is the thermal equilibrium 
number density (before freeze-out), < av > is thermally averaged value of the total annihilation 
cross section times the relative velocity of the annihilating neutralinos. The Hubble term takes care 
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Figure 2: The gaugino fraction of the lightest neutralino as function of mo and mi/2 for tan P=10 (left) and 
tan (3=50 (right). 
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Figure 3: Dominant tree diagrams for Dark Matter annihilation. Note that the amplitudes of the graphs shown 
at the top are proportional to the mass of the final state fermion, while the Higgs exchange is proportional to 
tan/3 for d-type quarks and _Z/tan/3 for up-type quarks. This implies that light fermion final states can be 
neglected and at large tan /3 the bottom final states dominate. 




Figure 4: The neutralino annihilation total cross section for i&n (3 =35 as function of the center of mass momenta 
in GeV of the neutralinos for quark -, lepton - and W'^W" final states, as calculated with CalcHEP ]61^ . Note 
the helicity suppression at low momenta for light fermions. 




Figure 5: The neutralino annihilation cross section for pseudoscalar Higgs exchange for bottom and top final 
states as function of tan /3 , as calculated with CalcHEP Ifilf . 
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Figure 6: The thermally averaged annihilation cross section times velocity for neutralino annihilation as 
function of mg and rni/2 for tan/3= 50 and bb, ti, W~^W~ , and tt final states (clockwise from top left). The 
results were calculated with DarkSusy JEBl- 
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Figure 7: The thermally averaged total annihilation cross section times velocity for neutralino annihilation as 
function o/ too and mi/2 for tan/3= 5 (left) and 50 (right).) The neutralino mass equals ~ 0.4toi/2 in the 
CMSSM, so in the plots the neutralino varies from ^0 to 400 GeV along the front axis. Note the strong decrease 
of the cross section for heavier SUSY mass scales and the different vertical scales due to the strong increase of 
the cross section with tan/3. 



5 



of the decrease in number density because of the expansion, while the first term on the right hand 
side represents the decrease due to annihilation and the second term represents the increase through 
creation by the inverse reactions. 

At temperatures below the mass of the neutralinos the number density drops exponentially. The 
annihilation rate T =< av > drops exponentially as well, and if it drops below the expansion 
rate, the neutralinos cease to annihilate. They fall out of equilibrium (freeze-out) at a temperature 
of about m^/25 [HT! and a relic cosmic abundance remains. 

For the case that < av > is energy independent, which is a good approximation in case there is 
no coannihilation, the present mass density in units of the critical density is given by 0: 

n = . (2) 

^ Pc <(JV> ' ^ ' 

One observes that the present relic density is inversely proportional to the annihilation cross section at 
the time of freeze out, a result independent of the neutralino mass (except for logarithmic corrections). 
For the present value of Q^^h? = 0.1 the thermally averaged total cross section at the freeze-out 
temperature of must have been 3 • 10"^"^ cm^s~^. This can be achieved only for restricted 

regions of parameter space in the MSSM, as will be discussed in the next section. 



3 Dark Matter Predictions from Supersymmetry 



Supersymmetry 1^2] presupposes a symmetry between fermions and bosons, which can be realized in 
nature only if one assumes each particle with spin j has a supersymmetric partner with spin \j — 1/2| 
(|j — 1/2| for the Higgs bosons). This leads to a doubling of the particle spectrum. Unfortunately the 
supersymmetric particles or "sparticles" have not been observed so far, so the sparticle masses must be 
above the limits set by searches at present accelerators. Obviously SUSY cannot be an exact symme- 
try of nature; or else the supersymmetric partners would have the same mass as the normal particles. 
The mSUGRA model, i.e. the Minimal Supersymmetric Standard Model (MSSM) with supergravity 
inspired breaking terms, is characterized by only 5 parameters: mo, mi/2, tan/3, sign(//), Aq. Here 
mo and mi/2 are the common masses for the gauginos and scalars at the GUT scale, which is deter- 
mined by the unification of the gauge couplings. Gauge unification is still possible with the precisely 
measured couplings at LEP 59 . The ratio of the vacuum expectation values of the neutral compo- 
nents of the two Higgs doublets in Supersymmetry is called tan /3 and ^o is the trilinear coupling at 
the GUT scale. We only consider the dominant trilinear couplings of the third generation of quarks 
and leptons and assume also Aq to be unified at the GUT scale. The constraints on the supersym- 
metric parameters space are practically independent of due to a coincidence from the constraints 
from the b — > rate and the lower limit on the Higgs mass of 114 GeV j^H]- The absolute value 

of the Higgs mixing parameter ^ is determined by electroweak symmetry breaking, while its sign is 
taken to be positive, as preferred by the anomalous magnetic moment of the muon |59j . 

The GUT scale masses are connected to low energy masses by the Renormalization Group Equa- 
tions (RGE), as shown in Fig. ^ The running masses of the gauginos at low energy obey the simple 
solutions of the RGE: 



ai{t) 

am 

Numerically at the weak scale {t = 2ln{MGUT / Mz) = 66) one finds (see fig. P): 



Miit) = ^rn,/2- (3) 



(4) 
(5) 
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The gluinos obtain corrections from the strong coupling constant 03; therefore they grow heavier 
than the gauginos of the SU{2)l U{1)y group. Since the Higgsinos and gauginos are all spin 1/2 
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particles and are equal in all other quantum nubers, the mass eigenstates are in general mixtures 
of them, which are called generically charginos (neutralinos) for the mixture of the super symmetric 
partners of the charged (neutral) gauge bosons and charged (neutral) Higgs bosons. The Majorana 
neutralino and Dirac chargino fields can be written as: 

/ ^ \ 

X- 

\ m I 

while the mass matrices can be written as j52j : 




(Ml -Mz cos 13 sinOw sin /5 sin 6*14/ \ 

M2 M z cos (3 cos 9w —Mzsm.j3cos9w 

—Mz cos j3 sin 9 w Mz cos (3 cos 9w — /Lt 

\ Mzsm.j3 sm.9w — sin/? cos — / 



(7) 



^ -\V2MwcosP /. ) 

The last matrix leads to two chargino eigenstates Xi2- "^^^ dependence on the parameters at the 
GUT scale can be estimated by substituting for M2 and /x their values at the weak scale: Mi{Mz) = 
2M2{Mz) ~ 0.4mi/2 and ii{Mz) ~ 0.63//(0). 

From Fig. ^ it can be seen that the mass parameters in the Higgs potential, mi and m2, are 
driven negative, largely because of the large Yukawa couplings of the third generation of quarks and 
leptons. This leads to radiative electroweak symmetry breaking (EWSB), so the Higgs mechanism 
in Supersymmetry needs not to be introduced ad hoc, as in the Standard Model, but is caused by 
radiative corrections. The running is only strong enough for top masses between 140 and 200 GeV and 
if the starting value y/j^ + m^ at the GUT scale is large enough, which in practice implies ^ > M2. 
From the mass matrices and El it is clear that for Mi^2 < A* the lightest chargino is wino-like with a 
mass given by M2 if the mixing is neglected and similarly the lightest neutralino is bino like with a 
mass given by Mi ~ O.5M2. 

In practice, there is some mixing and the neutralino mass eigenstates are linear combinations of 
the weak eigenstates, i.e. 

X° = iVil^ > +iV2|W'3 > +N3\H^ > +Ni\Hl > . 

The gaugino fraction A^^ + N2 is nevertheless close to one, especially if the diagonal elements are 
large compared with the off-diagonal elements proportional to Mz- This is demonstrated in Fig. |2 

The interaction of the sparticles with normal matter is governed by a new multiplicative quantum 
number called R-parity, which is needed in order to prevent baryon- and lepton number violation. In 
GUT theories quarks, leptons and Higgses are all contained in the same supermultiplet, which allows 
couplings between quarks and leptons. Such transitions, which could lead to rapid proton decay, 
are not observed in nature. Therefore, the SM particles are assigned a positive R-parity and the 
supersymmetric partners have a negative one, which can be related to the known conserved quantum 
numbers of spin S, baryon number B and lepton number L by = (_i)3-B+/+25 j^gq^ij-jng R-parity 
conservation implies that at each vertex one needs two supersymmetric particles, from which it follows 
that: 

• The rapid proton decays involving vertices with only one sparticle do not occur. 

• Sparticles can be produced only in pairs, e.g. pp qgX or e'^ + ^ /^^A • 
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• The heavier sparticles can decay to hghter ones, hke e — > 67 or g ^ qg, but the Lightest 
Supersymmetric Particle (LSP) is stable, since its decay into normal matter would change R- 
parity. 

• The LSP has to be neutral to be a good candidate for Dark Matter. 

• The interactions of particles and sparticles can be different. For example, the photon couples 
to electron-positron pairs, but the photino does not couple to electron-positron - or selectron- 
spositron pairs, since in these cases the R-parity would change from -1 to +1. 

• The LSP is weakly interacting with normal matter, since the final state has to contain the LSP 
again, so its interaction with quarks would only be elastic scattering by e.g. Z-, Higgs or squark 
exchange. 

Consequently, the LSP is an ideal candidate for Dark Matter, since it has all the properties of a 
Weakly Interacting Massive Particle (WIMP), namely it is neutral, heavy and weakly interacting, so 
it will form galactic haloes. 

The neutralinos can annihilate through the diagrams, shown in Fig. 13 The main features of the 
amplitudes have been indicated below the diagrams: 

• The annihilation into fermion-antifermion pairs is proportional to the fermion mass in the limit 
v^O, which is the important case in the present universe at a temperature of a few Kelvin. This 
can be understood as follows: The neutralino is a Majorana particle, so it is its own antiparticle. 
In addition it has spin 1/2, thus obeying Fermi statistics, which implies it cannot have identical 
quantum numbers. Furthermore at low velocity it annihilates into an s-wave state, which implies 
the spins have to be antiparallel (just like for electrons in the hydrogen s-wave). Therefore also 
the spins in the final state have to be antiparallel, which leads to an amplitude proportional to 
the fermion mass to account for the required helicity fiip^^. Consequently heavy final states 
are enhanced at low momenta, as demonstrated in Fig. ^ Note that at higher momenta not 
only s- waves contribute and the helicity suppression disappears, so during the time of freeze-out 
all final states were produced. Note that these arguments are only valid for the diagrams with 
sfermion and Z-exchange. For the Higgs exchange the proportionality to the final state fermion 
mass arises from the Yukawa coupling. These cross sections were calculated with the program 
package CalcHep|filj. 

• The second important point concerns the tan P dependence: the diagram via pseudoscalar 
Higgs exchange is proportional to tan/3 for down- type quarks and 1/tan/? for up- type quarks. 
This implies that at large tan /? (tan /? > 5) the b-quark final states are enhanced over t- 
quark final states, as shown in Fig. 13 The amplitudes of the Higgs exchange and the Z- 
exchange have an opposite sign. Since the top final states have a large amplitude for Z-exchange 
(amplitude oc mass), they are additionally suppressed by the negative interference with the t- 
channel amplitudes. The cross sections for various final states are shown as function of mo, Tni/2 
in Fig. El The strong increase of the total annihilation cross section as function of tan /3 is 
demonstrated in Fig. [3 in the mo,mi/2 plane, which also shows the strong dependence on 
the SUSY masses: only low SUSY scales and large tan/3 lead to cross sections of the order of 
magnitude required by the WMAP data (Eq. EJ. 

• As shown in Fig. |21 the amplitude for pseudoscalar Higgs exchange is proportional to iViA^3^4, 
i.e. it requires that the lightest neutralino has both bino- and Higgsino components, which 
implies that the diagonal elements in the mass matrix should not be too large compared with 
the off-diagonal elements proportional to Mz ■ So unless one tunes tan /3 and the SUSY masses 
such that one hits the resonance (m^ ~ 2m^), in which case very small Higgsino admixtures 
are enough, one needs relatively light neutralino masses. 
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Figure 8: Propagated abundances at 7.5 GeV/nucleon, as calculated by Galprop in comparison with data. From 
Ref. l^. 



In summary, the annihilation cross section becomes large for large tan (3 and is dominated for 
tan /3 > 5 by the s-channel pseudoscalar Higgs exchange into bb quark pairs. Fig. d shows that values 
of tan P around 50 yield the annihilation cross sections required by WMAP, given in Eq. |2 Regions 
of coannihilation at smaller tan /? are allowed by WMAP data, but have too small annihilation 
cross sections to explain the deficiencies in the positron, antiproton and gamma ray fluxes, as will 
be shown in the last section. Before discussing the global fits, the cosmic ray fluxes from nuclear 
interactions are discussed. 

4 Cosmic Rays generated by Nuclear Interactions 

The sources of charged and neutral cosmic rays are believed to be supernovae and their remnants, 
pulsars, stellar winds and binary svstems j53l I54j . Observations of X-ray and 7-ray emissions from 
these objects reveal the acceleration of charged particles near them. Particles accelerated near the 
sources propagate tens of millions of years in the interstellar medium where they can loose or gain 
energy and produce secondary particles and 7-rays. The spallation of primary nuclei into secondary 
nuclei gives rise to rare isotopes. Nuclear interactions produce not only matter, but also antimatter, 
like antiprotons and positrons. The latter originate mainly from the decay of charged pions and kaons. 

The detailed studies of cosmic rays teach us about the production and propagation in the universe. 
The gammas can deliver information over intergalactic distances, while the charged particles pro- 
pagate mainly on galactic distances. Secondary nuclei are produced in the galactic disc, from where 
they escape into the halo by diffusion and Galactic winds (convection). They may gain energy by 
"diffusive" reacceleration in the interstellar medium by the 2nd order Fermi acceleration mechanism, 
i.e. on average more interstellar gas clouds from opposite directions are hitting a given nuclei than 
"comoving" gas clouds. In elastic collisions this leads on average to an energy gain, thus depleting 
the low energy part of the source spectrum. Long-living radioactive secondaries tell how long they 
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Figure 9: The B/C ratio as an example of secondary /primary nucleon ratios for various Galprop models 
in comparison with data. The dashed lines are the model with diffusive reacceleration, the solid ones for 
diffusion plus convection, and the dotted ones for plain diffusion. The lower (upper) curves of each kind are 
the interstellar (solar modulated) ones. The diffusive reacceleration curve (dashed) provides the best fit to the 
B/ C ratio, but due to the large diffusion coefficient required, this leads to a deficiency in the flux of antiprotons 
(right). From Ref. Q/. 



survive in the halo before interacting in the disc again, thus determining the size of the halo. The 
gas density and acceleration time scale can be probed by the abundances of the K-capture isotopes, 
which would decay via electron K-capture in the interstellar gas. 

A global fit to all this information allows one to build a model of our galaxy. Analytical and 




Energy, MeV 

Figure 10: The gamma-ray energy spectrum of the inner Galaxy as calculated by Galprop ^ in comparison 
with EGRET data. Clearly, there is an access of data above 1 GeV. 
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semi-analytical models often fail when compared with all data. Therefore advanced models incorpo- 
rating nuclear reaction networks, cross sections for production of antiprotons, positrons, 7-rays and 
synchrotron radiation, energy losses, convection, diffusive reacceleration, distribution of sources, gas 
and radiation field etc. are needed. 

In addition, the distributions of matter and antimatter in the interstellar medium (ISM) are 
modified locally by affects of the solar activity and magnetic fields inside our solar system, e.g. from 
the planets. Gleeson & Axford modelled the periodically varying solar activity with a typical half 
cycle of 11 years^ by a radial solar wind in which the charged particles loose kinetic energy depending 
on their rigidity R and on the distance r from the sun. On this time scale the incoming flux from the 
Galaxy does not vary and the problem reduces to an adiabatical deceleration by the solar wind with 
a dependence only on the radial coordinate. This can be solved analytically: 



where J(r, t) is the measured differential flux at a distance r from the sun for particles with energy 
E and mass -Eq/c^, J(inf, £") is the incoming flux to the solar system and E? = (E'^ — EQ)/(Ze)'^ is 
the rigidity with R' = R{E'). The energy loss can be parametrized by the solar modulation parameter 
^{t) as E = E' — |Ze|$(t), where \Ze\ is the absolute charge of the particle; <I>(t) varies between 350 
and 1500 MeV depending on the solar cycle. The solar modulation shifts the particle spectrum to 
lower energies, but the effect is only noticeable for rigidities below 10 GV. Recent determinations of 
the local interstellar flux (LIS) from the modulated (=measured) ones for electrons, positrons and 
protons can be found in Ref. |64j . 

The most complete and publicly available code for the production and propagation of particles in 
our galaxy is the Galprop code |53l I54j . It provides a numerical solution to the transport equation 
including a cross section database with more than 2000 points, source functions, density distributions, 
etc. The cross section tables include all possible cross sections: p + p, p + He, p + N, He + N, N + N, 
where all nuclei up to the heaviest ones (Ni) are considered. Fig. |H1 shows the composition of 
the primary and secondary nuclei, as calculated by Galprop in comparison with data. Clearly, the 
production of secondary nuclei is well described. 

Fig. ini shows the spectrum of the Boron over Carbon (B/C) ratio, which shows a characteristic 
depletion at low energies. Since Boron is a purely secondary produced nuclei, while Carbon is primarily 
produced, the depletion at low energy is a sensitive handle on the question of diffusive reacceleration 
and solar modulation. As shown, the modulation effects the spectra mainly at kinetic energies below 
10 GeV/nucleon. In order to reproduce the sharp peak in the ratios of secondary to primary nuclei 
without any unphysical breaks in the energy dependence of the diffusion coefficients and/or the 
injection spectrum, the diffuse reacceleration with a rather large diffusion coefficient is needed 0. 
But this leads to too few antiprotons, as shown by the dashed line on the right hand side in Fig. El 
The possible way out of the discrepancy between the B/C ratio and too few antiprotons was suggested 
by Strong and Moskalenko: a "fresh" local unprocessed component at low energies of primary nuclei, 
thus decreasing the B/C ratio and allowing for a smaller diffusion coefficient. Also too few gammas are 
generated by Galprop, as shown in Fig. Il()| which would also need either a harder nucleon spectrum 
or a harder electron spectrum, but this would need spatial variations which make the spectrum in 
our local region unrepresentative of the large scale average |H1 • 

However, an alternative explanation may be the annihilation of neutralinos, which increases the 
yield of gammas, antiprotons, and high energy positrons, but does NOT effect the B/C ratio. This 
goes exactly in the direction of solving these discrepancies between Galprop and present data simul- 
taneously, as will be shown in the next section by a global fit to all data. 

cycle can vary between 8 and 14 yr. 




(9) 
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Figure 11: The halomodel used in this paper in comparison with other halomodels in the literature (left) and the 
integrated fraction of the fluxes as function of the distance from the center of the galaxy (right). One observes 
that more than 50% of the gamma rays originate from a region of less than 0.2 (1-5) kpc from the center for a 
detector subtending a solid angle of 10~^ (10~^) sr, while the antiprotons and positrons reach the 50% at 1.9 
and 4-1 kpc, respectively. 




Energy [GeV] Energy [GeV] 

Figure 12: The flux of antiprotons and positrons normalized to the gamma flux from neutralino annihilation 
before (left) and after (right) diffusion with the halomodel and diffusion parameters described in the text. 
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Figure 13: Boost factors for gamma rays, positrons and antiprotons (from left to right) as function of tuq and 
^^1/2 for ta.ii (3 = 5 1 . Note that the boost factors are similar for all three fluxes, which is a strong constraint for 
any model. 



5 Global Fits to positrons, antiprotons and gamma rays 

Trying to disentangle the contributions from nuclear interactions and neutralino annihilation to the 
antimatter fluxes and gamma rays is in practice not easy. Ideally one would like to implement the 
neutralino annihilation as a source function in the Galprop code, so the antimatter and gamma rays 
from nuclear interactions and annihilation would be transported through the galaxy in an identical 
way. However, this numerical code is too slow to be used in a fit program. Therefore, we used 
the second best possibility, namely using the publicly available code DarkSusy |55^ for neutralino 
annihilation, which has semi-analytical solutions to the diffusion equation and includes the important 
energy losses for positrons. We changed the diffusion parameters and code in DarkSusy in such a 
way, that the results resembled as closely as possible the Galprop results. The main difference in 
the diffusion parameters between DarkSusy and Galprop lies in the fact that Galprop uses diffusive 
reacceleration, while DarkSusy does not. The diffusive reacceleration is needed to fit the B/C ratio 
and results in an almost order of magnitude larger diffusion constant with a much smaller energy 
dependence. The diffusion parameters and energy losses in DarkSusy were changed as follows: 




Figure 14: The relic density (left) and boost factor for the gamma rays as function of mo and mi/2. Note 
that for large mass scales the relic density is reduced due to significant coannihilation cross sections, which are 
not operative in the present universe, thus causing large boost factors. For small mass scales the boost factor 
is practically proportional to the relic density, since both are inversely proportional to the annihilation cross 
section. 
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The size of the galaxy is determined by a cyHndrical box with half height h=4 kpc and radius 
r=30 kpc, as in Galprop. The solar system is at a distance of 8.5 kpc from the center. Outside 
this cylinder the diffusion stops and the particles escape. 

The energy dependence of the diffusion constant was taken to be the same as in Galprop and 
given by: 

D = (3D^- {p/por, 

where P = v/c is the velocity, Dq = 42 • 10^^'' crv? j s is the diffusion coefficient at momentum 
Po = 1 GeV and q = 0.33. This diffusion coefficient is at least a factor 7 larger than the default 
one in DarkSusy and has a much weaker energy dependence, which is of importance mainly for 
antiprotons. 

The large diffusion constant implies that the galactic center starts to contribute significantly 
to the antiproton flux in the solar system. Therefore the effective thickness of the disk was 
increased from 0.1 to 0.5 kpc and the density of hydrogen in the disc was increased from 1 to 3 
GeV/cm~'^, thus increasing the interactions between antiprotons and protons in the disc. The 
density of hydrogen in the halo was neglected. 

The energy loss dE/dt = (^E"^ of the positrons (mainly by synchrotron radiation and inverse 
Compton scattering on star light) was doubled using C = 1.52 • 10^^ yr-^ GeV~^, following 
Galprop [HSl and [HH] . 

The Dark Matter halo profile is usually assumed to be of the Navarro, Prenk & White (NFW) 
type |67j . as supported by numerical simulations of galaxy formation |68l I69j . Effectively we 
have chosen for the Dark Matter density distribution an isothermal spherical symmetric profile: 



,'r 



p{r) = po-{-r' 1 + i-r 



a 



a 



7-/3 



where a is a scale radius and the slopes a, /3 and 7 can be thought of as the radial dependence 
at r ~ a, r >> a and r << a, respectively. At large distances we expect a 1/r^ dependence for 
a flat rotation curve, while at small distances more a 1/r dependence is needed. For definiteness 
we use the default (a,/3, 7) =(1,3,1) for a scale a = 10 kpc, but e.g. (1.5,2,1) for a=2 kpc yields 
practically identical results. At a distance less than 10~^ kpc from the center the profile is kept 
constant to avoid any singularity. The density po is adjusted such that the local halo density is 
in the range 0.2-0.8 GeV/cm^, as required by the rotation curve of our galaxy j51j . The halo 
density increases quite steeply towards the center, so most annihilations will take place in the 
center of the galaxy, thus producing there gamma rays, positrons and antiprotons. Different 
halo profiles are compared in Fig. 1111 The percentage of gamma rays coming from the center 
depends on the halo profile and on the solid angle, as shown on the right hand side of Fig. 
111! The positrons only reach the detector from a much smaller region due to the higher energy 
losses, while the antiprotons come from an intermediate region, as shown in Fig. ^2 as well. 
These curves were obtained by calculating the flux for a halo profile truncated after a given 
radius. 

The gamma rays can travel over large distances without loosing energy, but they will arrive 
at the detector only, if they were emitted along the line of sight. Antiprotons and positrons 
on the other hand change direction during the propagation along the magnetic field lines and 
by collisions, so they can arrive at the detector even if they were not emitted along the line of 
sight. This causes a larger acceptance for the antiprotons and positrons in comparison with the 
gamma rays, as shown in Fig. 1121 Here the gamma flux at the Earth - per unit of time, 
surface and solid angle - from neutralino annihilation in the galactic halo can be written as: 

- ' Plds = —^^, <J>, (10) 



dt dS dO. 4tt m| Ji ^.s. 47r m| 



14 



where is the number of photons per annihilation (about 38 for b-quark final states) and 
< J > is the averaged value of the integral of the neutralino density squared along the line of 
sight. Since the density peaks at the center, the averaged value depends on the solid angle over 
which one averages. For the chosen profile < J > is about 20000 for a solid angle of 10~^ sr 
towards the center of the galaxy and 1000 for dO, = 0.17 sr, which is the solid angle used for 
the EGRET data. 

• Diffusive reacceleration is an essential ingredient of Galprop, of special importance for the an- 
tiprotons. DarkSusy has neither the antiproton background spectrum from nuclear interactions 
nor the possibility of reacceleration. Therefore we used the parametrization of the antiproton 
spectrum from Simon, Molnar and Rosier |70j . The effect of the reacceleration is a flattening of 
the antiproton spectrum, which was simulated in DarkSusy by a larger value of the solar mod- 
ulation constant <I>. With the parametrization from Simon et al., the large diffusion coefficient 
and ^ = 1200 MeV we could reproduce the Galprop spectra for data from the years 1997 and 
1998 sufficiently well. We will use only data from this part of the solar cycle. 

The annihilation cross sections in DarkSusy differs considerably from the one in Micromegas |HH] 
for large tan /3. This was traced back to a different width of the pseudoscalar Higgs boson. Micromegas 
follows the usual procedure of minimizing the Higgs potential at a scale given by ^yfnf^rn^, which 
minimizes the higher order corrections. Furthermore, the pseudoscalar Higgs width is determined by 
the Yukawa couplings, of which the most important one in our case is the one given by m;,. However, 
the QCD corrections to the width lead to expressions of the type mb{fJ.){l + cias{n) /it ln(//^/m^) + 
the large logarithms are minimized by choosing for the renormalization scale = rriA- After using 
the running b-mass in DarkSusy to calculate the width, good agreement in the cross sections was 
found for the large tan /? region, where pseudoscalar Higgs exchange dominates. 

The relic density was calculated with Micromegas as well, since this program is particularly 
suited for large tan /3, where the Higgs exchange and its loop corrected width are important and 
it furthermore incorporates all coannihilation channels. Micromegas was interfaced to the program 
package Suspect, which calculates the low energy SUSY masses from the GUT scale parameters |56j . 
The problem of calculating the relic density at large tan /3 [Jlj is caused by the fact that both the 
pseudoscalar Higgs mass and its width have large corrections: the mass at tree level is determined 
by ml + evaluated at the electroweak scale, but as can be seen from Fig. ^ the running of 
becomes steep at low energies. Since m"^ becomes negative for large tan /3, the tree level mass becomes 
small or even negative and a positive final mass is obtained by radiative corrections. Therefore the 
different calculations give quite a spread in masses ft2j and depend sensitively on the scale, where 
the Higgs potential is minimized. 

The following data were used in the fit: 

• Gamma ray data from the galactic center in the angular range 330° < i < 30° and — 5° < 6 < 5° 
from the EGRET space telescope, which has been taking data for about 9 years on the NASA 
Compton Gamma Ray Observatory (CGRO). We use the data as presented in Ref. j73j . 

• Positron data from AMS 174) and HEAT 

• Antiproton data from BESS in the years 1997 and 1998 |76| 

In order to see if the deficiencies in the galactic data on positron, antiprotons and gamma rays 
can be fitted by the contributions from neutralino annihilation, the following strategy was persued. 
A minimization was performed between the combined data and the sum of the annihilation signal 
and background. The normalizations of the fluxes of positrons, antiprotons and gamma rays from 
neutralino annihilation, the so-called boost factors, were taken as free parameters in the fit for a 
given neutralino mass, since the annihilation rate is proportional to the Dark Matter density squared, 
so any clustering of dark matter can increase the absolute normalization. The normalization of the 
background was not varied, since it is well determined by other data, which Galprop used to fit 
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the diffusion parameters. The Galprop data describes well the positrons and gamma rays in the 
region where the neutralino annihilation does not contribute, but fails for the high energy parts of the 
spectra. The antiprotons are too low over the whole spectrum for the diffusion parameters determined 
from the B/C fraction, as mentioned in the previous section. 

The function for gamma rays was defined as Y^iUi^i - Tf/{f^aif + (1 - f^if/al, where the 
sum runs over all data bins with errors CTj and T is the sum of the parametrized background from 
nuclear interactions, as calculated by Galprop and the contribution from neutralino annihilation, as 
calculated by DarkSusy with the modifications mentioned above. /-^ is a common normalization factor 
for all data points with a systematic error given by cr„. Normally was set to one, but sometimes 
it was left free in the fit in order to study the effect of possible correlations between the data points. 
Similar functions were defined for the antiprotons and positrons and the total was simply the 
sum of the contributions for gamma rays, antiprotons and positrons, since no correlations exist. 

The x^ function was minimized by the Minuit package [77| and the fit was repeated for all values 
of the SUSY parameters mQ,mi/2,tan (3, where the mass scales were varied between 200 and 1000 
GeV and tan /? was varied between 50 and 55. The sign of the parameter was chosen positive, as 
preferred by the data on Br(b Xsj) and the anomalous magnetic moment of the muon §^ — 2 |59j . 

The boost factors for positrons, antiprotons and gamma rays are shown as function of mo and 
mi/2 for tan/3=51 in Fig. IIISI They were all close to each other over most of the parameter space, 
which is absolutely non-trivial, if one considers the large enhancements for positrons and antiprotons 
from the diffusion as compared to gamma rays, as shown before in Fig. 1121 This is a strong indication 
that the deficiencies in the prediction of the antiprotons, positrons and gamma rays are indeed due 
to the missing dark matter annihilation with a standard dark matter halo profile. One expects a 
similar boost, since the fluxes all originate from quite a large region of the galaxy, as shown before 
m Fig. ini so the averaging over the dumpiness should be similar for all. The boost factors are 
determined by the present annihilation cross section at a temperature of a few Kelvin^ while the relic 
density is determined by the cross section during freeze out, which is about ra^j^h or 10^^ K. At 
the high temperatures coannihilations can become important, especially at high mass scales, so the 
annihilation cross sections can become large, but the boost factors can be large as well due to the 
small annihilation cross section in the present universe. This is demonstrated by a comparison of the 
relic density and the boost factor in Fig. 1141 where at large mg and the relic density is small 
due to coannihilation, but the boost factors are large. Small boost factors are only obtained at small 
mass scales for which the neutralino annihilation cross sections are large (see Fig. [7|. 

The regions of parameter space allowed by the WMAP data are plotted in Fig. El for different 
values of tan /3. It is clear that for tan /3 ~ 50 only a small region is allowed, if in addition the 
electroweak constraints from the Higgs mass and h ^s7, and the requirement that the LSP is 
a neutral particle have to be fulfilled. Values of tan /3 below 50 are excluded completely, if one 
wants to be consistent with all constraints and if one requires in addition that the boost factors 
are below 10. The last requirement implies that dumpiness can enhance the annihilation signal 
by at most a factor of 10, which is the value suggested by simulations of galaxy formation and it 
excludes the regions of coannihilations between staus and neutralinos. The boost factors are strongly 
correlated with the value of the local halo density p-^. To obtain a conservative excluded region by 
the requirement that the boost factor is below 10, was set to its maximum allowed value of 0.2- 
0.8 GeV/cm^. Coannihilations occur if the staus are nearly degenerate in mass with the neutralinos, 
which happens next to the region labelled "excl. LSP" in fig. 1151 In this region the stau is the Lightest 
Supersymmetric Particle (LSP) and next to it the stau is almost degenerate with the neutralino, so it 
cannot decay into a neutralino and tau. In this case the stau is practically stable and can annihilate 
with a neutralino into a tau plus photon. This coannihilation reduces the relic density to values 
required by the WMAP data. For values of tan /? above 50 there are quite a range of neutralino 
masses allowed, since for larger values of tan /3 one hits the resonance of pseudoscalar Higgs exchange, 
in which case much heavier neutralinos are allowed, as shown in Fig. (bottom left). 

The results of two good fits with two different neutralino masses (207 and 427 GeV) are shown 
in Figs. 1161 1171 and 1181 The parameters are summarized in Table ^ Fig. shows that these 
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parameters are close to masses at the lower and upper range of possible masses. As indicated in 
the figures the boost factors for antiprotons, positrons and gamma rays are all similar for the NRW 
halo profile discussed before and the improves significantly with the inclusion of Dark Matter in 
the fits. The x^/d.o.f. is reduced from 110/38 for the "background-only" fit to 34/35 for the fit 
including neutralino annihilation. This corresponds to an increase in probability from about 10^^ to 
0.5, as shown in Table [21 together with the results for the individual spectra. For the antiprotons the 
increase in probability is the least significant, as expected, since the shape of background and signal 
are similar. This causes also a strong correlation between the normalization of the data and the boost 
factor. For the gamma rays and positrons this correlation is less severe because of the different shapes 
of signal and background. Leaving all the normalizations free in the fit does not change the results 
significantly for normalization errors below 10 %, since mainly the antiprotons are affected, but they 
anyhow do not have a strong statistical significance. They are nevertheless important for the analysis, 
since they check the expectation of similar boost factors for all species and constrain the analysis by 
the requirement of not overproducing the antiprotons. 

It is interesting to see what are the prospects to observe the neutralinos in other channels or by 
direct detection for the points of parameter space selected by the present analysis. For a neutralino 
mass of 207(427) GeV the flux of muons with energy above 1 GeV from neutralino annihilation inside 
the sun are only 8 (0.3) km~'^ yr'^, while the flux from the earth is even smaller for the SUSY and 
halo parameters considered in this analysis. These fluxes were calculated with DarkSusy [55) . 

The cross section for direct Dark Matter detection was calculated with DarkSusy [55j as well. For 
the heavier masses the cross section rapidly decreases, but the projected sensitivity of future direct 
detection experiments may be sufficient, if the neutralino mass happens to be at the lower range 
allowed by the present analysis, as demonstrated in Fig. 1191 

Unfortunately, the high cross section deduced from the 6.3a solar modulation signal observed in 
the DAMA experiment [ZH| seems difficult to reconcile with the present analysis, since the diagrams for 
neutralino-nucleon scattering are by crossing symmetry similar to the ones for neutralino annihilation, 
so such a large cross section in direct detection would produce an excessive amount of positrons, 
antiprotons and gamma rays. 
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Figure 15: The region of relic density allowed by the WMAP data. The upper row is for ianj3=51 and = 
(left) and A^j = mo (right), which shows that the role of Higgs constraint (dotted line) and b X^j constraint 
(solid line) are interchanged for the different values of the trilinear coupling, but the lower limit on TOi/2 is not 
very sensitive to Aq. The light shaded (blue) area is the region allowed by WMAP and the contours of larger 
Q,h^ are indicated by the dashed lines in steps of 0.05. The second (third) row show the same information for a 
larger region for tan (3—51, 52 (I. and r.) ( 53 and 55 (I. and r.)). The excluded regions, where the stau would 
be the LSP or EWSB fails or the boost factors are above 10 are indicated by the dots. The black dots indicate 
the resonance region, where |m^ — 2mj(-(,| < 10 GeV. Jbrtan/3 > 52 the acceptable values for the relic densities 
are for TOi/2 values above the resonance region. 
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Table 1: mSUGRA parameters with the corresponding supersymmetric particle spectrum and some predicted 
observables for an LSP mass of 207 and 427 GeV. The cross sections were calculated by DarkSusy ]55}/ . which 
agree with the cross sections from Micromegas ]58l/ for these parameters after the modifications to DarkSusy 
mentioned in the text. The particle spectrum was calculated by Suspect ^5.6^ and Feynhiggsfast \57Jl . 
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59.4/18 
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0.48 


14.3/17 


0.59 


Antiprotons 


15.9/13 


2.6-10-1 


6.6/12 


0.89 


6.7/12 


0.88 


Gammas 


34.5/7 


1.4-10-^ 


9.5/6 


0.15 


12.3/6 


0.06 
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110/38 


6.7-10-y 


33.7/35 


0.53 


33.3/35 


0.55 



Table 2: The and probabilities of the fits for fluxes from background only, i.e. from nuclear interactions 
(labelled "b" ) and for fluxes including "signal" contributions from neutralino annihilation (labelled "b+s" ) for 
two neutralino masses. 
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Figure 16: Gamma ray spectrum with contributions from nuclear interactions (grey /yellow) and neutralino 
annihilation (dark/red) for a neutralino mass of 207 (left) and ^21 GeV (right). 
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Figure 17: Antiproton spectrum with contributions from nuclear interactions (grey/yellow) and neutralino 
annihilation (dark/red) for a neutralino mass of 207 (left) and ^27 GeV (right). 
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Figure 18: Positron spectrum with contributions from nuclear interactions (grey/yellow) and neutralino anni- 
hilation (dark/red) for a neutralino mass of 207 (left) and 4^7 GeV (right). 




WIMP Mass [GeV] WIMP Mass [GeV] 

Figure 19: The spin independent (left) and spin dependent (right) cross section for direct WIMP detection 
on protons in comparison with existing (dashed lines) and projected (solid lines) limits and the cross section 
expected for neutralino masses of 144 '"^'^ ^^^5 GeV (squares). A neutralino mass as low as 144 GeV is at 
the border of compatibility with the constraints from the Higgs mass and b Xg^, so a mass between the two 
masses indicated is more likely. 
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6 Conclusion 



The galactic models can very well describe the nuclear reactions, especially the ratio of secondary to 
primary produced nuclei and the small concentration of long-lived radioactive nuclei. These obser- 
vations fix the diffusion parameters, especially the need for diffusive reacceleration as an additional 
term in the diffusion equation was emphasized by Moskalenko and Strong. However, the diffusion 
parameters from such a global view produce too few antiprotons, too few hard gammas and to a lesser 
extent to few hard positrons. Up to now many investigations on the fluxes of positrons, antiprotons 
and gamma rays have been performed, either to study the contributions from neutralino annihilation 
or to study "Local Bubble" type contributions to defeat the deficiencies of the standard propaga- 
tion models. However, it is shown that Dark Matter annihilation within the Constrained Minimal 
Supersymmetric Model can solve these deficiencies simultaneously. 

This is the first time that these fiuxes are studied in a global fit using the best available diffusion 
parameters without artificial constructs to remedy the deficiencies of the galactic models and add to 
the contributions from nuclear interactions the ones from neutralino annihilation in the Constrained 
Minimal Supersymmetric Model (CMSSM) with all the known constraints from electroweak precision 
data and the relic density from the WMAP satellite. Actually, all we need from the supersymmetric 
model to fill up the deficiencies is a neutralino with a thermally averaged annihilation cross section 
times relative neutralino velocities of about 10~^^ err? j s for the present universe, i.e. without coan- 
nihilation. This cross section can be smaller, if the halo distribution of the neutralinos is not smooth, 
but clumpy, thus enhancing the annihilation rate by a certain "boost factor", which simulations of 
galaxy formation show to be of "the order of a few". The boost factors for antiprotons, positrons 
and gamma rays were left as independent free parameters in the fit, but they all come out to be 
similar for a given halo profile and given neutralino mass. This result is non-trivial given the very 
different mean free paths of positrons, antiprotons and gamma rays and provides a strong support 
for neutralino annihilation as a common origin of the deficiencies. The needed large cross sections 
are obtained in the CMSSM only for relatively light neutralinos in the mass range of 150 - 400 GeV 
and tan /? > 50. At large tan /3 the annihilation is dominated by pseudoscalar Higgs exchange into hh 
quark pairs. The upper neutralino mass is limited by the rather arbitrary requirement that the boost 
factors are below 10 and the lower mass limit is given by the electroweak constraints at large tan/3. 
The probability of the global fit to galactic fiuxes of positrons, antiprotons and gamma rays increases 
from 10~* to 0.5 by including the contribution from neutralino annihilation, which can be interpreted 
in case of gaussian errors as an improvement by about 6 standard deviations. These facts, statistical 
significance of the global fit combined with similar boost factors for positrons, antiprotons and gamma 
rays, provide strong experimental evidence for the supersymmetric nature of Dark Matter. 
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